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Compressor  Modeling  for  Engine  Control  and  Maintenance 

Carl  Beckey,1  Roy  Hartfield,2  Mark  Carpenter3 
Auburn  University,  Auburn,  AL 

Abstract 

This  work  describes  a  method  for  dynamically  updating  the  compressor  map  for 
an  advanced  turbofan  engine.  The  implementation  of  the  method  is  illustrated  using  data 
from  a  transport  aircraft  powered  by  four  turbofan  engines.  The  compressor  map  is 
generated  through  projection  pursuit  regression  and  non-linear  approximation,  using 
Bernstein  polynomial  basis  expansion  regression.  This  approach  allows  for  prediction 
options  which  would  be  useful  to  a  controller  algorithm.  Results  include  the  compressor 
maps  and  recommendations  for  implementing  the  proposed  approach  in  an  engine 
management  system. 

Introduction 

The  refined  turbofan  engine  offers  unmatched  performance  and  efficiency  for 
heavy  lift  transport  applications.  This  paper  focuses  on  further  refinements  to  the 
operation  of  the  compressor.  As  air  is  ingested  by  the  compressor,  it  travels  through  a 
complex  array  of  rotating  and  stationary  airfoils,  pressure  bleed  valves  and  ducting.  To 
analyze  this  flow  and  its  effects  on  engine  performance,  real-world  data  has  been 
captured  from  a  turbofan  engine  and  a  parametric  model  of  the  operating  conditions  has 
been  developed.  The  data  captured  during  a  range  of  operations  has  been  used  to  develop 
a  compressor  map.  Purely  mathematically  surrogate  approaches  to  modeling  compressor 
performance  have  been  developed1,  however,  in  this  case  a  parameterized  modeling 
approach  is  used  because  of  the  higher  likelihood  of  capturing  some  of  the  physics  of  the 
compressor  performance  with  the  model.  This  data  has  been  reduced  using  Bernstein 
polynomials  to  generate  a  set  of  speed  lines  and  a  complete  response  surface  for  the 
compressor.  This  approach  can  be  implemented  dynamically  to  provide  a  continuously 
updated  performance  map  available  to  the  engine  controller  across  the  entire  operating 
range  for  the  engine.  As  components  wear  and  deteriorate  throughout  the  life  of  the 
engine,  operating  conditions  can  be  adapted  by  the  dynamic  updating  of  the  compressor 
map  response  surface.  This  provides  full  utilization  of  the  operating  envelope  from  the 
time  an  engine  is  placed  on  wing  until  it  is  removed,  allowing  for  maximum  efficiency, 
performance,  and  safety. 
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Experimental  Procedures 


Data  captured  from  each  of  the  four  engines  of  a  heavy  transport  aircraft  were 
used  to  conduct  this  analysis.  This  family  of  engines  has  a  6:1  bypass  ratio  and  a 
maximum  overall  pressure  ratio  of  just  over  30.  As  seen  in  Figure  1,  compression  is 
provided  by  a  dual  spool  compressor.  The  low  pressure  compressor  (LPC)  consists  of  a 
fan  and  four  compressor  stages,  while  the  high  pressure  compressor  (HPC)  consists  of  a 
set  of  variable  pitch  inlet  guide  vanes  (IGVs)  and  12  compressor  stages,  of  which  the  first 
four  have  variable  pitch  stators.  Bleed  air  is  pulled  off  four  stages  throughout  the 
compressor  to  support  a  number  of  systems  including  aircraft  pressurization,  air- 
conditioning,  anti-ice,  inert  gas  generation,  active  clearance  control  and  engine  stability. 
Data  is  captured  at  a  number  of  locations  throughout  the  engine  and  combined  with 
general  aircraft  operating  conditions  to  provide  detailed  performance  feedback.  In  the 
development  of  the  analytical  model,  total  pressure  in  pounds  per  square  inch  from 
stations  2,  3  and  4.9,  total  temperature  in  degrees  Celsius  from  stations  2  and  4.9,  weight 
of  fuel  flow  in  pounds  per  hour  and  RPM  as  a  percent  of  maximum  for  the  HPC  and  LPC 
have  been  utilized.  The  data  has  been  recorded  from  multiple  aircraft  over  a  period  of 
several  days,  during  which  time  over  ten  thousand  data  points  were  captured.  The  data 
points  include  aircraft  operation  at  a  number  of  different  airfields  and  conditions  both  in 
the  United  States  and  abroad.  Furthermore,  the  data  points  cover  a  range  of  operating 
conditions,  including  climb,  cruise,  descent,  landing,  air-to-air  refueling,  low-level  flight, 
high  speed  descents  and  multiple  different  airdrop  profiles. 


Figure  1:  Turbofan  station  numbering  diagram. 
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Data  Reduction 


The  measurements  as  outlined  above  have  been  utilized  to  generate  a  compressor 
map  and  response  surface  for  the  compressor  operation.  In  order  to  do  so,  it  was 
necessary  to  first  calculate  the  mass  flow  rate  through  the  core  of  the  engine.  For  this 
calculation,  the  conditions  just  downstream  of  the  turbine  were  utilized.  At  this  location, 
the  total  pressure  was  cross-referenced  to  the  inlet  conditions  to  check  for  choked  flow 
conditions,  as  such  the  local  velocity  was  assumed  to  be  sonic.  As  a  result,  all  data  points 
where  the  pressure  ratio  was  less  than  unity  were  rejected.  Furthermore,  to  improve 
response  accuracy,  all  data  at  and  below  ground  idle  RPM,  were  rejected.  The  remaining 
data  was  further  reduced  as  follows.  The  total  temperature  aft  of  the  turbine  was  utilized 
to  calculate  the  local  ratio  of  specific  heats,  which  allows  for  the  calculation  of  the  mass 
flow  parameter.  For  this  calculation,  the  specific  heat  at  constant  pressure  was  first  found 
using  Equation  1 } 

Cp  =  0.2269807 e0-000097247*7!  (1) 

The  ratio  of  specific  heats  is  then  given  by  Equation  2. 


The  mass  flow  parameter  is  then  calculated  for  the  point  just  downstream  of  the  turbine 
where  the  flow  is  assumed  to  be  sonic,  Mach  equal  to  one.  Given  the  ratio  of  specific 
heats  and  the  Mach  number,  the  mass  flow  parameter  is  given  by  Equation  3. 3 


MFP 


Thyjri 
A  Pt 


M  1  + 


Zzl  luzVa-r) 


The  hydraulic  diameter  can  then  be  utilized  to  find  the  area  at  this  station.  Substituting 
the  mass  flow  parameter  into  Equation  4,  in  combination  with  the  area,  total  pressure  and 
total  temperature  at  this  station,  the  mass  flow  rate  aft  of  the  turbine  can  be  calculated. 

rh  =  *  MFP  (4) 

Vrt 

The  mass  flow  rate  through  the  compressor  core  is  given  by  the  mass  flow  rate  at  the 
turbine  exit  minus  the  mass  flow  rate  of  the  fuel,  converted  from  pounds  per  hour  to 
pounds  per  second.  In  order  to  isolate  the  performance  of  the  compressor,  the  data  must 
be  reduced  such  that  the  affects  of  atmospheric  temperature  and  pressure  are  normalized. 
This  is  accomplished  by  non-dimensionalizing  the  data,  as  shown  in  the  following 
equations.  Corrected  mass  flow  rate  is  given  by  Equation  8. 
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The  pressure  ratio  for  the  compressor  is  simply  given  by  Equation  9. 


ratio 


Pt2 


(7) 

(8) 


(9) 


Utilizing  Equations  7  through  9,  the  compressor  map  was  generated  by  equally 


distributing  an  array  of  RPM  values  across  the  operating  range  of  the  collected  data. 
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Figure  2:  Scatter  plot  of  corrected  data  organized  by  RPM  groups. 


Non-linear  Regression  and  Linear  Basis  Expansion 

Following  the  notation  in  Hastie,5  et  al.  (2008),  let  X  6  HP  denote  a  real  valued  random 
input  vector  (mass  flow  rate)  and  Y  EH  a  real  valued  random  output  variable  (i.e., 
pressure  ratio),  with  joint  distribution  Pr(X,  7).  We  seek  a  function  f(X)  for  predicting  Y 
given  values  of  the  input  X.  Based  on  a  loss  function  L(X,f(X))  for  penalizing  errors  in 
prediction, 
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L(Y,f(X))  =  (Y-f<iX))2 


we  would  like  to  find  the  function  that  minimizes  the  expected  prediction  error  (EPE) 
defined  as  EPE(f)  =  EL(Y,  /(X)).  For  squared  error  loss,  it  is  easy  to  see  that 

EPE(f)  =  EL(Y,f(X )) 

=  E(Y-f(X ))2 

=  j[y-  f(x)]2  Pr (dx,  dy ) 

=  ExEy]x([Y  -  f(X)]2 \X) 

For  fixed  values  of  X,  we  have 

EPE(f)  =  Ey{x([Y  ~  f(X)]2  \X) 

And  we  see  that  the  EPE  is  minimized,  for  each  x,  at 

fix)  =  E(Y\X  =  x). 

Note  that  f(x)  is  known  as  the  regression  function  and  can  be  linear  or  non-linear.  If  the 
function  is  non-linear,  we  may  still  be  able  to  estimate  it  using  a  linear  expression 
through  linear  basis  expansion  functions.  In  this  approach,  we  expand  the  design  space 
using  linear  basis  expansion,  such  as  is  done  in  polynomial  regression.  Denote  hjm(X )  : 
MP  —>  M  by  the  mth  transformation  of  X,  m  =  1, .  .  .  ,  M.  We  model 

p  p  M  j 

f(X )  =  £/>(*;)  =  2  Z  Pjmhl '»■(*')■ 

7  =  1  7  =  1  772=1 

which  represents  linear  basis  functions  in  X.  Note  that  the  expression  is  linear  in  the 
parameters,  /?;m,  m  =  1, ... ,  Mj  and  j  =  1, ... ,  p.  Some  examples  are 

a.  hm(X)  =  Xm,  m  =  1, ... ,  p,  which  is  the  basis  for  multiple  linear  regression. 

b.  hm (X)  =  X2  or  hm(X)  =  XjXk,  which  allows  for  a  higher-order  Taylor 
expansion  and  represents  the  basis  for  quadratic  response  surface  regression  models. 

c.  hm(X)  =  XlfXlf  ...  XpP,  ij  =  0,1  ,j  =  1,  ...,p,  which  allows  for  higher-order 
polynomials  and  interactions.  If  we  let,  h^X)  =  Xj,j  =  1,2, ... ,  p,  h2(X)  = 

Xj  Xj,  i  <  j  ,h2  (X)  =  Xt  XjXk,  i  <  j  <  k, ...  ,hp (x)  =  XxX2  ...  Xp,  we  will  have 
2P  —  1  terms  generated.  In  the  case  of  p=3,  we  have,  if  we  include  an  intercept 
term, 


fix)  =  Po  +  PlXl  +  p2x2  +  P3X3  +  Pax1x2  +  p5xlx3  +  Psx2x3  +  p7xlx2x3 
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In  this  paper,  we  follow  the  above  regression  approach  using  the  Bernstein  polynomials 
as  the  basis  expansion  functions.  These  are  described  below. 

The  data  was  then  grouped  by  RPM  and  a  curve  was  fit  to  each  RPM  group.  A 
number  of  different  fits  were  explored  in  order  to  generate  the  best  response  for  the  data 
and  included  a  linear  fit,  a  least-squares  polynomial  fit  and  a  Bernstein  polynomial  fit. 
To  preserve  the  characteristics  of  the  compressor  operations  while  allowing  for  a 
dynamic  fit  to  the  data,  a  Bernstein  polynomial  was  selected.  This  approach  allows  the 
function  to  remain  fixed  at  the  predicted  surge  and  choke  lines  while  allowing  the 
shaping  characteristics  of  the  polynomial  to  accurately  represent  the  intermediate 
operating  conditions  to  the  order  n.  The  Bernstein  polynomial  typically  defined  over  the 
interval  0-1  has  been  scaled  to  operate  on  a  continuous  function  over  the  interval  a  to  b, 
which  is  a  function  of  the  particular  RPM  being  evaluated.  For  the  Bernstein 
polynomials  defined  in  Equation  10  by 

<10) 

there  are  n+1  nth-degree  Bernstein  polynomials,  each  term  of  which  can  be  most  easily 
represented  by  the  Pascal  triangle  shown  in  Figure  3. 


Figure  3:  Pascal  triangle  of  individual  Bernstein  polynomial  terms  up  to  n=4. 

Before  the  Bernstein  polynomials  could  be  applied  to  the  data  set  it  was  first  necessary  to 
define  a  series  of  control  points.  The  control  points  for  a  number  of  speed  lines  are 
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Normalized  Pressure  Ratio 


represented  by  the  vertices  of  the  linear  spline  can  be  seen  in  Figure  4  and  Figure  5. 

Compressor  Scatter  Plot  for  Tail#  1 
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Figure  4:  Control  points  and  resulting  speed  line  curves. 
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The  control  points  were  found  by  dividing  each  speed  line  into  equally  spaced  segments 
and  are  shown  by  the  pink  and  dark  blue  points  in  Figure  5.  This  data  is  further  refined 
by  selecting  only  the  25%  which  most  closely  match  the  RPM  being  analyzed.  These 
data  points  are  further  separated  into  lower  and  upper  values,  pink  and  dark  blue 
respectively,  from  which  a  lower  and  upper  value  is  computed  as  the  mean  of  their 
respective  data  points,  depicted  by  the  red  in  Figure  5.  Finally  the  control  point  is 
linearly  interpolated  from  the  resultant  upper  and  lower  values  for  the  given  interval. 
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Figure  5:  RPM  segment  for  the  computation  of  a  single  control  point. 
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The  complete  set  of  control  points,  Pi9  for  each  speed  line  is  used  to  compute  the  Bezier 
Curve,  P(t)9  using  Equation  11. 

P(i)  =  J%=0PlBiin(t)  (11) 

The  complete  set  of  control  points  and  the  resulting  curves  can  be  seen  in  Figure  4  and 
the  complete  compressor  map  is  shown  in  Figure  6. 
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Compressor  Map  for  Tail#  1 
Engine#  3 
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Figure  6:  Single  engine  complete  compressor  map. 

The  Bernstein  polynomial  curves  along  with  the  normalized  data  from  Equations  7,  8  and 
9  were  then  further  utilized  to  interpolate  between  RPM  groups  using  a  cubic  spline  and 
build  a  complete  response  surface  for  the  compressor. 
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Figure  7:  Single  engine  complete  response  surface. 
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Results  and  Discussion 


The  operating  data  from  the  engines  of  a  single  aircraft  were  processed  as 
discussed  above  and  the  results  are  summarized  in  the  following  figures.  In  Figure  8,  a 
number  of  the  features  of  a  compressor  map  become  evident  as  the  corrected  values  are 
plotted.  First,  the  speed  lines  begin  to  appear  as  the  various  colors  representing  different 
RPMs  group  together  and  emerge  from  the  data  set.  Second,  overall  performance  of  the 
compressor  and  engine  are  easily  discernable.  Third,  a  line  connecting  the  data  points 
along  the  far  left  hand  side  of  the  graph  approximates  the  surge  line  for  the  compressor. 
Operations  above  the  surge  line  result  in  a  mass  flow  rate  which  is  insufficient  to 
maintain  the  given  pressure  ratio.  The  bleed  valves  at  stages  5,  14  and  17,  along  with  the 
variable  IGVs  and  stators  within  the  engine,  are  used  to  relieve  the  pressure  and  prevent 
cavitations  which  result  in  sudden  loss  of  thrust  and  possibly  even  engine  damage. 
Fourth,  a  line  connecting  the  data  points  along  the  far  right  hand  side  of  the  graph 
approximates  the  choke  line  where  the  compressor  is  reaching  its  mass  flow  limit. 
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Figure  8:  Compressor  map  scatter  plots  with  RPM  grouping,  approx  surge  and 

choke  lines  for  installed  engines. 


Figure  8  shows  a  large  variance  of  pressure  ratio  and  mass  flow  rate  for  each 
depicted  interval  of  HPC  RPM.  This  is  most  likely  a  result  of  the  numerous  different 
configurations  the  compressor  bleed  valves  and  stator  angles  can  have,  giving  rise  to 
differing  flow  parameters.  In  order  to  account  for  these  variations  and  develop  a  more 
accurate  parametric  model  of  the  compressor  operations,  further  data  pertaining  to  bleed 
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settings  is  required.  Figure  9  depicts  each  speed  line  showing  a  good  approximation  of 
the  operating  conditions  due  to  the  endpoints  being  constrained  by  the  surge  and  choke 
lines  and  the  shaping  terms  of  the  Bernstein  polynomial  fitting  the  values  in  between. 


Compressor  Map  for  Tail#  2  with  Bernstein  Polynomial  Speed  Lines  of  Degree  50 
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Figure  9:  Speed  lines  fit  using  a  Bernstein  polynomial  for  installed  engines. 

Bernstein  polynomials  were  applied  to  the  speed  lines  and  a  complete  response 
surface  has  been  generated.  The  results  can  be  seen  in  Figure  10.  The  response  surface 
provides  an  analytical  method  to  compute  pressure  ratio,  given  any  combination  of  HPC 
RPM  and  mass  flow  rate.  The  response  surface  has  been  validated  by  calculating  the 
coefficient  of  determination  (R2).  The  residual  for  each  of  the  points  has  been  calculated 
and  the  resulting  sum  of  squares  of  residual  (SSE)  along  with  the  total  sum  of  squares 
(SST)  and  R2.  The  results  for  each  surface  are  tabulated  versus  the  data  from  its  own  and 
the  accompanying  three  engines  from  the  aircraft  in  Table  1.  Furthermore,  a  visual 
representation  of  the  goodness-of-fit  is  provided  by  Figure  11,  which  gives  two  different 
perspectives  of  the  response  surface  overlaid  with  the  captured  data  points  for  its  own 
engine.  Both  the  graphical  depiction  and  the  calculated  values  for  the  goodness-of-fit 
suggest  this  is  a  well  behaved  model  which  accounts  for  almost  99%  of  the  variance  in 
pressure  ratio  capable  of  accurately  modeling  performance  of  this  family  of  engines 
operating  under  a  given  set  of  conditions. 
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Figure  10:  Compressor  with  Bernstein  polynomial  response  surface,  installed  engines. 


Table  1:  Goodness-of-fit  results  for  each  response  surface  vs.  each  engines’  data  capture. 
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Figure  11:  Response  surface  and  captured  data  overlay. 

Further  performance  analysis  is  conducted  by  comparing  the  response  surfaces  of  each 
engine  against  each  other  engine  on  an  aircraft.  For  this  analysis,  a  jet  with  engines  of 
significantly  different  operating  life  times  is  selected.  Figure  12  shows  the  margin  plots 
for  this  aircraft  and  clearly  shows  engine  #2  having  a  much  lower  margin,  approximately 
-3°C,  than  any  other  engine  on  this  aircraft. 
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Figure  12:  Engine  operating  temperature  margin  as  a  function  of  time. 

When  each  point  on  the  response  surface  is  subtracted  from  the  matching  point  on  each 
of  the  other  response  surfaces,  the  following  figures  are  generated.  An  overhead 
perspective  of  the  resulting  plots  clearly  illustrates  portions  of  the  operating  regime  where 
compressor  ratios  are  higher  or  lower  than  the  comparative  engine.  Figure  14  is  of 
particular  interest  as  it  shows  the  results  of  a  low  margin  engine.  Here,  it  can  be  seen  that 
the  compressor  is  operating  with  higher  pressure  ratios  under  most  conditions  and,  in 
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particular,  along  the  surge  line  when  compared  to  the  others  as  a  result  of  its  increased 
operating  demands  in  order  to  produce  the  same  exit  pressure  ratio  as  the  other  engines. 
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Figure  13  Engine  #1  compressor  performance 
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Figure  14:  Engine  #2  compressor  performance. 


14 

Approved  for  public  release;  distribution  is  unlimited. 


AFFTC-PA-1 1253 


Normalized  Mass  Flow  Rate  Normalized  Mass  Flow  Rate  2  Normalized  Pressure  Ratio  Normalized  Mass  Flow  Rate 


Compressor  Difference  Surface  for  Tail#  3 


Engine#  3  minus  Engine#  1 


N2  RPM 


Engine#  3 


N2  RPM 


xIO 


0.6  0.9 
Drmalized  Mass  Flow  Rate 


1 0.9 

LL 

CO 

c§  0.8 


Engine#  3  minus  Engine#  2 


/ 


1 


2 

1 


0 


n  0.7 


o  0.6 


1  1.05  1.1 

N2  RPM 


x  10 


Engine#  3  minus  Engine#  4 


Figure  15  Engine  #3  compressor  performance. 
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Figure  16  Engine  #4  compressor  performance. 


15 

Approved  for  public  release;  distribution  is  unlimited. 


AFFTC-PA-1 1253 


Conclusion 


Data  from  a  modem  compressor  has  been  captured  from  a  wide  range  of  missions 
and  operating  conditions,  providing  a  complete  platform  upon  which  to  conduct  an 
analysis.  The  data  was  first  normalized  to  eliminate  variation  from  atmospheric 
conditions,  and  then  it  was  used  to  calculate  the  corrected  mass  flow  rate  and  the  overall 
pressure  ratio  of  the  compressor.  The  resultant  data  was  further  reduced  to  generate 
speed  lines  using  a  Bernstein  polynomial  approach  and  a  compressor  map  was  generated. 
The  compressor  map  was  then  further  developed  to  provide  a  complete  response  surface 
for  the  compressor.  The  surface  provides  a  function  by  which  pressure  ratio  can  be 
derived  from  HPC  RPM  and  compressor  mass  flow  rate.  The  response  surfaces  provide 
the  surge  and  choke  lines  of  the  compressor.  The  accuracy  of  the  model  has  been 
determined  by  inspection  of  the  data  and  calculation  of  the  coefficient  of  determination 
which  suggests  that  the  model  accounts  for  almost  99%  of  the  variation  in  the  response 
variable,  pressure  ratio.  Further  analysis  of  the  response  surface  can  be  used  to  provide 
optimal  and  safe  control  and  operation  over  the  life  of  the  engine. 
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